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I .  Introduction 

The  synthetic  aperture  radar  on  SeaSat  Satellite  has  provided  high  resolu¬ 
tion  synoptic  images  of  the  earth's  surface.  However,  the  same  geographical 
regions  viewed  at  different  look  angles  sometimes  provided  different  images. 

This  especially  was  the  case  of  the  images  in  the  urban  areas  surrounding  the 
Los  Angeles  areas  (1).  In  these  images  the  areas  of  the  bright  returns  did  not 
correspond  while  certain  areas  of  low  backscatter  were  common  to  both  the  . 
images.  A  study  of  the  causes  of  this  variation  in  the  imagery  is  the  subject 
matter  of  the  present  report. 

A  possible  model  for  the  study  of  backscatter  from  urban  areas  due  to  the 
synthetic  aperture  radar  is  to  consider  it  as  a  periodic  gridded  structure  with 
the  streets  and  the  buildings  forming  the  grid.  The  average  dimensions  of  the 
grid  for  our  study  must  have  dimensions  large  compared  to  the  wavelength. 

As  part  of  this  task  a  literature  survey  was  conducted  and  the  principal 
relevant  references  were  found  to  be  due  to  DeSanto  (2),  Jordan  &  Lang  (2), 
Holford  (4),  Uretsky  (5),  and  McCammon  and  McDaniel  (7).  Uretsky's  method  which 
is  based  on  the  Helmholtz  equation  reduces  the  scattering  from  a  periodic 
surface  to  an  infinite  set  of  linear  equations  of  the  form  [A]  [X]  =  [Y]. 

DeSanto  also  arrived  at  the  same  result  using  a  different  approach.  For  solving 
the  infinite  set  of  equations,  the  method  of  reduction  was  employed.  In  the 
method  of  reduction  a  truncation  of  the  infinite  series  is  employed  and  the 
correctness  of  the  solution  will  depend  on  the  order  of  matrix.  However,  for 
equations  of  the  first  kind,  there  does  not  exist,  at  present,  a  criteria  for 
truncation  which  will  assure  the  convergence  of  the  solution. 

Hoi  ford's  method,  like  that  of  Uretsky,  is  based  on  the  Helmholtz  integral 
formula  and  results  in  a  Fredholm  integral  equation  of  the  second  kind  for  the 
scattered  Intensity.  This  reduces  his  integral  equation  to  the  form 
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[I]  +  [A] 


[X] 


[V] 


and  rigorous  proof  exists  for  the  convergence  of  the  solutions.  We,  therefore, 
plan  to  use  the  approach  due  to  Hoi  ford  to  study  the  scattering  from  a  periodic 
surface.  The  only  numerical  results  available  using  Holford's  approach  is  due 
to  McCammon  and  McDaniel  and  here  they  consider  only  relatively  shallow  struc¬ 
tures,  viz,  for  Kh  1,  where  h  is  the  height  of  the  structure,  K  *  2tt/a  and  A 
is  the  periodicity. 

Section  II  of  the  report  contains  a  detailed  review  of  the  relevant  liter¬ 
ature.  Section  III  has  the  results  of  the  backscatter  calculations  for  a 
typical  geometry  suited  for  our  observations  using  the  exact  method  due  to 
Hoi  ford.  Section  IV  has  the  conclusions  and  the  recommendatl ons  for  future  work 
on  the  problem. 


1.  SeaSat  Views  of  North  America,  The  Caribbean,  and  Western  Europe  with  Imaging 
Radar,  JPL  Publication,  60-67,  California  Institute  of  Technology,  Pasadena,  CA, 

2.  J.A.  DeSanto  (1981),  Rad.  Sci.  16,  p.  1315. 

3.  A.K.  Jordan  and  R.A.  Lang  (1979),  NRL  Report  8284. 

4.  R.L.  Hoi  ford  (1981)  JASA  70  p.  1116. 

5.  J.L.  Uretsky  (1966),  Ann.  Phys.  33.  p.  400. 

6.  D.F.  McCammon  and  S.T.  McDaniel  (1985),  JASA  78  p.  149. 


II.  Review  of  the  Literature 


The  general  methods  of  treating  periodic  structures  can  be  broadly  classi¬ 
fied  into  two  general  approaches.  The  first  approach  are  the  so-called  approxi¬ 
mate  methods  and  the  second  are  the  more  exact  methods.  Sections  1 1 -1  and  1 1 -2 
review  briefly  the  approximate  methods,  while  Sections  I I -3  through  I I -5  review 
the  more  exact  methods.  In  Section  I I -6,  some  of  the  details,  that  have  to  be 
considered,  for  the  numerical  implementation  of  the  calculations  for  the  Hoi  ford's 
method  are  discussed. 


I I -1  Kirchoff  Method 

The  general  Kirchoff  Method  is  based  on  Helmholtz  integral  equation. 
Consider  a  plane  wave  incident  on  a  rough  surface.  The  scattered  field  at  any 
observation  point  P  is 


(1) 


where  the  integral  is  over  the  surface  and  E  is  the  field  on  the  surface. 
4>  is  the  free  space  Green's  function  (exp( i kR^ )/R" ) , 


where  R'  is  the  distance  from  the  origin  to  P  and  k  *  (2ir/x).  In  the  above 
equation,  E  is  an  unknown.  The  Kirchoff  approximation  assumes 

E  *  (1  +  R)E1  and  (f£)$  =  (1  -  R)E1 (k-S)  (3) 

where  n  is  the  normal  to  this  surface  at  P  and  R  is  the  reflection  coefficient 
of  a  smooth  plane.  In  the  case  of  an  one  dimensional  surface  where  L  =  c(x), 
and  R  *  +1 ,  Beckman  (1963)  arrives  at  the  following  for  the  scattering  coef- 
ficient,  p _  ’ _ 

A  Beckmann  and  A.  Spizzichino  (1963).  "The  Scattering  of  Electromagnetic  Waves 
from  Rough  Surfaces",  Pergammon  Press. 
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p  ( ©  1  * ©2 ^  3  (^/Eq)  s 


’l+cos(01+0g)' 

1 

coso^+cosag 

A 

a/2 

-A/2 


eir*vdx 


(4) 


0^  is  the  direction  of  the  incident  field,  Eg  is  the  field  reflected  in  the 
direction  62  and 

r*7  =  ~  [(sine^  -  sin0g)x  -  (cose^  +  cos0g)c(x)]  (5) 

This  approximation  is  valid  when  A>>\  and  the  radius  of  curvature  of  the  ir¬ 
regularities  is  large  compared  to  the  wave  length.  Consequently,  this  method 
will  break  down  when  there  are  sharp  edges. 


I 1-2  Rayleigh  Method 

In  the  Rayleigh  Method,  the  field  is  expressed  as  a -sum  of  plane  waves 
whose  directions  correspond  to  the  scattering  angles  and  for  a  periodic  surface, 
the  scattering  direction  will  be  given  by  the  grating  equation.  So  for  a 
periodic  surface  c(x) 

c(x)  ■  hcosKx  (6) 


The  grating  equation  for  the  field  is 


sine2m  *  sin0g  +  m  3  °»  i1 »  +2,... 


(7) 


The  total  field  in  the  region  z>c(x),  E,  (eq.  8),  is  the  sum  of  the  incident 
field,  E-| ,  and  the  reflected  field,  Eg, 


E  •  E,  *  E2 


.e<Vf+  *  A 

...  m 

m»-« 


(8) 


where 

*2m  *  k(ex  s1ne2m  +  V080*’ 

A  A 

<1  =  k(ex  sin0Q  -  ez  coseQ) 

AAA 

r  *  (xex  +  yey  +  zez) 

and  e  ,  e  and  e  are  the  unit  vectors  in  the  x,y  and  z  direction.  For  the  case 
x  y  z 

of  the  Dirichlet  boundary,  E|z_?(x)  =  0.  So,  from  Equations  (7)  and  (8),  we 
have 

e-ik;(x)cose0  a  _  j  A  e+imkx  +  ikc(x)cose2m 
m=-®  m 

Expanding  both  sides  of  Equation  (9)  in  terms  of  Bessel  functions,  Beckman  et 
al .  (1963)  arrives  at  the  expression  for  the  coefficient  Am 

J„(kKcos90)  *  j  (-<)”  \  Wkhcose2n,>  (10) 

m=-® 

From  this  set  of  equations,  the  coefficients  Am  are  determined,  which  are  used 
in  Equation  8  to  obtain  the  field. 

Uretsky  (1966)  notes  that  the  major  error  in  the  solution  given  in  Equation 
8  is  Its  implication  that  it  describes  the  field  everywhere  above  the  surface. 
This  has  also  been  pointed  out  by  Lipman  (1953)  that  the  solution  does  not  take 
into  account  both  up  and  down  going  waves  excited  in  the  "valley"  of  periodic 
surface.  We  will  now  describe  two  methods  of  solving  the  problem,  both  based  on 
the  integral  equation  formulation,  starting  from  the  Helmholtz  equation,  which 
will  take  into  account  both  the  upgoing  and  downgoing  waves  in  the  periodic 


structure. 
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I I -3  Uretsky's  Method 

This  section  describes  the  method  proposed  by  Uretsky  ("1966)  for  solving 
the  boundary  value  problem.  For  simplicity  we  consider  a  two  dimensional 
problem  and  the  scalar  wave  equation.  The  problem  reduces  to  the  solution 
of  the  two  dimensional  equation 


3  r(x,z)  .  a^P(x.z) 


2 

k> 


3X 


3z‘ 


01) 


subject  to  boundary  conditions  which  can  be  either  soft,  i.e.,  P(x,z)  =  0  on  the 

boundary,  or  hard,  i.e.,  *  0,  or  an  impedence  boundary  condition,  i.e. 

3  P 

—  +  AP  *  0  on  c(x),  where  A  is  a  constant.  Specifically,  we  will  consider  a 
surface 


z  ■  c(x)  *  hcosKx  where  K  »  (2tt/a)  (12) 

The  field  is  now  the  sum  of  the  incident  field  Pfnc(x,z)  and  the  scattered  field, 
P,  given  by  the  Helmholtz  integral. 


P(X,Z)  *  Pinc(x,z)  + 


P(r") 


3Hn(1)(k|?-F'|) 


3n' 


-H^tklr-r'l)  -^p-]ds  (13) 

where  r'  *  (x',  c(x'))  and  the  integral  is  taken  along  the  boundary.  In  the 
case  of  a  soft  boundary,  (P(r')  ■  0) 


where 


P(r)  *  P1nc(r) 


i|»(x')  H^(k|r-r"l)  dx' 


(14) 


^(x) 


I  if.  ds  I 

k  an  dxj 


z*c(x) 


£♦<-<*> 


if.] 

3XJz-q(x) 


(15) 


As  Equation  11  is  unaffected  by  a  shift  of  2ir/A  and  the  boundary  condition 
is  also  unaffected  because  it  has  a  periodicity  equal  to  2tt/a,  we  can  make  the 
following  expansion  for  ^(x). 

*(x)  =  l  *„e1kV  (16) 

n=-« 

where 

sine„  =  a  =  [a.  +  n  K/k) 
n  n  o 

cosen  ‘  Tn  ’  (1-“n2)V2;  “n-1 

=  K-1%2)1/21  °n>l  (17) 

Substitution  of  the  above  in  Equation  (14)  yields  (Holford  (1981)) 


where 


P(r)  =  Pjnc(r)  *  l  *„  I  PikV 


n=-<*>  m*-06 


i  Jktz,, 

■m ,n(z)  =  STT  j  Cm-n(kT) 


Cn(t)  .  1  j  e-itt(k)-1"l<xdx 
0 

Holford  (1981)  has  shown  that  the  integral  Im>n(z)  is  9iven  by 

W2>  ■  Cm-n<k^e<kV  z>t(x) 

Also  the  field  P(x,z)  can  be  written  as 

\  n  t  \  ,  r  q  ika„X  +  iky  z 

P(x,z)  =  Pinc(x,z)  +  Z  Rne  n 

n=-°° 

Inserting  Equation  21  into  Equation  18  and  comparing  with  Equation  22,  we 


obtain 


oo 

**m  ”  2y_  ^  ^n  ^m-n^Ym^ 


'm  n=-« 

We  now  proceed  to  determine  the  In  the  Equation  (18),  we  let  the  field 


8 


point,  r  approach  the  boundary. 

09 

P(r)|  ,  "  Pinc(r)!  ,  l  [  ^x>)  Ho1)(k|r-r>l)dx' 
z-c(x)  inc  z*c (x)  *  J  0 

—09 

If  we  now  apply  the  boundary  condition  that  P(r)  *  0  we  have 

09 

|  j  *(x')  H<1}  (k|r-r'|)dx'  -  -  P1nc(x) 

-o» 

Substituting  Equation  (16)  in  Equation  (25) 

09 

f  f  n'11  (k*>d*'  *  -p1nc(x) 

*  n 

*09 

where  y  «  [(x-x-)2  +  (c(x)  -  c(x'))2]^2 
Applying  the  operator 


(24) 


(25) 


(26) 


A 

H 


e-,kVdx 


to  both  sides  of  Equation  (26),  we  obtain  on  substituting,  p^nc(x,z) 
from  Eq.  8. 


(27) 


1k«  x-ikv  z 
o  o 


i  ••• 


where  V, 


m,n 


2lkVdx 


\  thJ11  (ky)e+ikotnx'  dx' 


(28) 


(29) 


This  set  of  linear  equations  (Equation  28)  Is  used  to  solve  for  4<n  which  is 
then  substituted  in  (Equation  23)  to  obtain  Rp. 

Note  that  the  set  of  linear  equations  Is  of  the  form  AX  *  Y.  The  procedure 
recommended  by  Uretsky  (1966)  is  to  truncate  the  set  at  some  value  n*N  and  then 
solve  the  linear  equation  by  matrix  Inversion.  However  as  Hoi  ford  (1981)  has 
pointed  out  there  is  no  guarantee  that  the  solution  for  AX  «  Y  will  converge  as 


Fig.  1.  Plane  wave  scattering  from  a  surface  s(x)  of  period 
L  and  depth  d.  The  closed  contour  (123451)  is  used  for  Green's 
theorem.  The  incidence  angle  is  0,,  and  8„  is  one  of  the  discrete 
scattering  angles.  Diagram  lines  are  slightly  offset  for  clarity. 


WlliiiiWl 
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iih-  or  will  converge  to  a  correct  solution.  This  method  is  valid  only  for 
equations  of  the  second  kind  represented  by  (I  +  A)X  *  Y  where  I  is  an  identity 
operator.  Holford's  method  of  overcoming  this  will  be  discussed  later  in 
Section  I 1-5. 

II-4  Method  of  Desanto 

Desanto  (1981)  has  proposed  another  approach  which  is  similar  to  that  of 
Uretsky.  For  a  sinusoidal  surface  shown  in  Figure  1  and  a  plane  wave  incident 
on  the  surface. 

*,(*,!)•  »1kV*1kV  (30) 

where 

9Q  ■  Incident  angle 

°o  *  sin0o» 

Yo  ‘  coseo 

The  reflected  field  is  then 


where 

sine,,  *  a„  ■  a„  +  n  x/A  ■  a,  +  n  K/k  n  «  0,  +1,  +2,... 
n  n  o  o  —  ~ 

cose,  » 
n  n 

The  Green's  function 

6±<x,z)  •  (ei(1ltoV-1koV))  (32) 

satisfies  the  wave  equation  in  the  region  z>o.  Me  now  consider  the  contour 
12345  in  Figure  1  and  apply  the  Green's  theorem  to  the  contour 

j  <♦*"-  61  UMu  •  0  (33) 

c 

and 


*{x,z)  =  4»fnc(x,z)  +  l  Anelko°nx  +  ikV  ( 

n=-“ 

If  we  now  let  23  and  51  shrink  to  zero,  then  the  integral  along  the  path  12 
becomes 

+A/2 

*  f  ti1*lcoYmel  ko^ao"am^X  -  iko‘Ym^Ane^k°^0tn"am^  +  l  koYoe<k°^0to?Qtn,^> 

-A/2  n 

“  ikJA*YneikoK+am)x}  dx  ( 

0“  n  n 
n 

For  G+  along  1-2,  Eq.  35  reduces  to 

‘l2  =  -4’lYin(A/x>smo  < 

and  for  G"  along  1-2,  Eq.  35  simplifies  to 


Hence 


*ia  * (4’^  r>  \ 


A  f-O  f°r  G 

mo. 

A_  for  G 


The  Integral  along  23  and  51  are  zero  since  we  have  allowed  2  to  approach  3  and 
5  to  approach  1.  The  equation  simplifies  to 


\  <♦  f- £  $  d"  *  0 

12  345 


1 


|®T  .  St  |t)  dy 
an  an'  M 


(^Ym  f”) 


for  + 


0  -Am  for- 


We  now  need  to  evaluate  the  left  hand  side  of  the  equation.  Let  us  assume 


Oirichlet  boundary  condition. 


-  |  «£  f£)  du  *  -(4irlYm  f~) 


“V  for  + 


o  A_  for  - 


Following  the  method  of  Uretsky,  we  now  replace  —  i|>(x,c(x))du  by  $(x)dx,  where  <j>(x) 
is  a  periodic  function  with  a  period  A. 

*(x)  -  |2-e1k0VB(x)  (42) 


where 


m 

B(x)  •  2  exp[-1krtdv„0  +  cos(2ir/A)x)]  l  inan  exp(in2irx/A) .  (42a 

o  o  „  ^  n 

n*-* 


Substituting  Eq.  42  to  the  right  hand  side  of  Eq.  41,  we  obtain 

|  S^(x.t(x))  du 
345 

A/2 

«  j  G±<x,c(x);  j^-e1koaoxB(x)dx 
-A/2 

Now,  let  c(x)  «  -  |<1  +  cos  and  let  y  ■  ~  x, 

|  G~(x ,c (x) )  1^-  du 
345 

■  V(1VV2)  f  flno+'m***'2  *  1n,y)C(y)dy 
* 


where 


C(y>  »  B (J-y) 

Using  equation  41,  we  obtain  for 


f  C (y)e(-1my  '  Wos*/2)  -  ,(1k0d/2ir.) 

J  ^ n 


mm® 
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and  for  G" 


f  C(y)  exp[-imy  +  i  ( k0dYfT,  cos  y/2)]  dy  3  -Am(4TriVV  exP^_ik0dV2^ 

'  -  TT 


(46) 


From  Eqs.  44(a)  and  (42a) 


C(y)  =  -  +  cosy)/2)  (1  +  l  iV"yl 


(47) 


o  n^o 

Substituting  for  C(y)  from  Equation  47  in  Equation  (45),  then  yields  the  set  of 
linear  equations  for  m^O 


J_(t)  +  l  a  J  ( t  )  3  0 
m  m  n  m-n  m 

Similarly,  substitution  in  Equation  (46)  gives 


(48) 


where 


Am  3  “  (-“)  imexp(iT^)  {J_(t^>  +  l  an  J  (t*)} 
m  y_  m  m  m  ‘  n  m-n  m 

*•  ■.*> 


(49) 


*  M 

i  *  (to  i  V  -§- 


We,  therefore,  solve  the  first  set  of  linear  equations  to  obtain  a„  which  is 

n 

then  used  in  Equation  49  to  obtain  the  A  . 

m 

Jordan  and  Lang  (1979)  solved  the  Equations  48  and  49  numerically.  They 

state  that  no  numerical  instability  was  experienced  and  the  stability  of  the 

2  Yn 

method  was  checked  by  testing  conservation  of  energy  (i.e)  X I Rn I  — =  1  for 

n  Yo 

real  orders  of  n.  Wirgin  (1980)  reports  that  the  method  of  Desanto  was  also 
used  by  Whitman  and  Schwerlng  (1977)  who  reported  no  instability  for  Kh  <_  1  and 
Kh  <_  2.36.  However,  in  their  experience  this  method  does  not  work  well  for 
Kh>ir .  Further  Wirgin  states  that  the  solution  will  become  unstable  for  Kh>u. 

A.  W(rg1n,  (1980),  JASA  68,  p.  692. 

G.  Whitman  and  S.  Schwering,  (1977),  IEEE.  PGAP  25,  p.  869. 


■ui 


y ■AC-'.  vv 
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Though  Jordan  and  Lang  reported  no  instability  in  their  numerical  computation, 
as  the  equations  are  of  the  form  AX  »  Y,  convergence  to  the  correct  solution 
cannot  be  guaranteed.. 


1 1-5  Method  Due  to  Hoi  ford 

Holford  (1981)  has  developed  a  method  to  overcome  this  problem.  This 

method  is  similar  to  the  approach  of  Uretsky,  where  we  start  with  the  two 

dimensional  boundary  value  problem  and  seek  a  solution  to  the  differential 

equation  satisfying  the  given  boundary  conditions  on  the  surface  and  the  radia- 

2  2 

tive  condition  at  x  +2  -*■  -. 

I 1-5-1 :  Surface  Satisfying  the  Oirichlet  Boundary  Condition: 

Applying  the  Helmholtz  Integral  formula  as  in  Section  I I -3  we  obtain  for 
the  Oirichlet  boundary  condition. 


OB 

p(r)  *  pinc(r)  +  |  j*(x')H0^(k|r-r'|)  dx' 


where 


*(*>  -  (r  £  §>.„ 


kk  an  dx'z»?(x) 

Let  the  normal  derivative  operator  ik  [-  yy  +  c'(x)  “J  be  applied  to  both 
sides  and  taking  on  the  limit  when  z  -*■  t(x),  we  obtain 


(50) 

(51) 


<»(x)  *  <|/{x)  +  y 


i'(x')  [-  |j+  c'(x)  |y]  H^(k|r-r'|)dx' 


(52) 


It  has  been  shown  by  Meecham  (1956)  that  the  Integrand  on  the  right  hand  side 
has  a  singularity  at  r  *  r'  and  this  integrates  to  a  term  ^(x)  in  the  limit. 


Thus,  we  have 

h>M  =  *(x)  -  4  ♦(X-)  J - 

2  4  |r-r'| 


[(c(x')  -  ? ( x) )  -  (x'-x)c'(x)]dx 


where 


♦(xW+1k*Vt-  fj  +  t'(x)  |j]e1kV  '  <kV 


■  -  tY0  +  a0C'(x)]eikV  •  1kV(X> 

We  therefore  obtain 

ik  7  nP^(k|r-r"| ) 

*(x)  +  ~2  U(x-)  - C(;(x')  -  c(x))  -c(x--x)c-(x)]dx- 

i  |r-r'| 

*  2*(x) 

As  explained  in  Section  II-III,  ip(x)  is  periodic  with  the  period  A.  So 


*(X>  *  f  *„e1kV  ■  l  v'k(<.0«<K/k)x 


■  e1kaoxIV,nKx 
n 

Substituting  for  ^(x)  in  the  above  expression 


I<.ne1kV  e1k  nx'(H(’)(k|?.?-|)/|r-r-|) 

n  n  ; 


CU(x')  -  c(x))  -  (x'-x)c'(x)]dx'  *  2i|»(x) 
We  now  apply  the  operator 

1  f  e‘ikV  dx 

•  t 


to  both  sides  of  the  equation  and  obtain 


where 


fit  +  Jib  U  *  2w 

vm  £vn  m,n  mi 


1  -ika  x  .  k 
—  e  m  dx  s- 

A  2 


Ika  X'  ^’(klf-F'l) 

e  n  - 


CU(x')  -  t(x))  -  (x'-x)t'(x)]dx" 

A 

1  [;<x)e-ikv 

o 

*-r|  <*0  w'(*>>e1k”°x  ‘ 1l<v(x)  d* 


This  set  of  linear  equations  (Eq.  59)  Is  now  solved  to  determine  <|>n  and  then  ij>(x) 
This  Is  then  used  to  obtain  P(r).  The  specular  reflection  coefficients  for 
various  orders  of  scattering  is  obtained  from  Equation  (23)  which  is  reproduced 


below 


"  57“  £  ^n  Sn-n^Ym^ 


I 1-5-2:  Surface  Satisfying  the  Neumann  Boundary  Condition: 

Applying  the  boundary  condition  |jj-  *  0,  in  Eq.  13,  we  obtain 

o» 

P(r)->  P1nc(r)  *  jf  |  «(x')  -3f  [H0(,)(k|f-?'|)](z,  ,  5{x-})  f^dx- 


where  $(x')  *  P(x',c(x')).  We  now  let  r  -+  [x,?(x)],  then  using  the  same 
method  as  was  used  in  deriving  Eq.  52,  we  obtain 


i 


\  «(x)  a  ♦(x)  + 


f  ♦(«') 


r-r 


r(c(x)  -  t(xO)  -  (x-x')c'(x')]  dx' 

U  7  H,(1)(M 

$(x)  =  2<j>(x)  +  -g-  <fr(xO  - - - C(x,x')  dx' 
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(64) 

(64a) 


where 


C(x,x')  =  [(c(x)  -  c(x'))  -  (x-x')e'(x')] 
u  »  [x-x')2  +  (x(x)-x(x1)2]1/2 


(65) 


where  <j>(x)  *  p^nc(x.  c(x)).  We  now  write 

♦(*)  *  I  *„e1k“"X 

n»-« 


Substituting  in  equation  (64a),  we  obtain 

oo 

£<t>ne*kanx  -  C(x,x')(£<j>ne1Kanx  ) 


ika^x'x  ik 


dx'  =  2<p(x) 


(66) 


(67) 


A 

Applying  the  operator  ^  j  e"ikctmxdx  to  both  sides,  we  obtain 


-i ka_x  . 
e  m  dx 


(I«ne1k“nx'}  if  -]  C(x,x') 


dx- 


r  - 


2  — 
4  A 


$(x)e”lkamx  dx 


** 

<b  -  T  d>  u  =  <b 
ym  -  n  m,n 


where  U„  „  is  given  by 
m,n 


U 


m,n  A  J 


1  f  e-ikV  dx 


ikHl0)(kw)  ika  X' 

I4-J— - e“lkV  C(x,x')  dx' 


(68) 

(69) 


(70) 
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The  set  of  linear  equations  obtained  can  now  be  solved  to  yield  $  which  then 
can  be  substituted  in  the  Equation  63  to  yield  P(r).  The  reflection  coefficient 
for  this  case  is  given  by  the  following  expression  (Holford). 


■j  (m-n)Ka_ 

R  3  s —  J  <p  C  (Kv  )(y  +  — r - ) 

m  Zy  L-  n  m-n  m  m  ky_ 


We  note  that  in  this  case,  the  equation  for  $n  is  in  form  (I  +  A)X  *  Y  and 
therefore,  convergence  is  guaranteed. 


I I -6  Numerical  Implementation  of  Holford  Method 


We  will  now  describe  the  numerical  implementation  of  the  method  due  to 
Holford,  as  we  will  be  discussing  the  results  In  the  following  Section  III. 
Specifically  we  will  consider  the  case  where  the  surface  is  a  sinusoidal  surface 
and  satisfies  Neumann  conditions.  We,  therefore,  have 

c(x)  *  h  cos  Kx 


The  linear  set  of  equations  to  be  solved  is 


<b  ■  I  4  V  3  4 
¥m  £  ¥n  m,n 


m  *  0,  +  1 ,  ... 


where  „  is  given  by 


A 

-  1 
i.n  a  _ 


j  -Ika  x  ik  f  ikax'  H1  (klr-r'U 
dx  e  m  — *-  e  n  - 

L  I?-?'!  z « dx) 


(t(x)  -  c(x')  -  (x  -  x')t'(x'))  dx' 


Where 


I r-r'|  «  [(x  -  x')2  +  (c(x)  -  r(x'))2] 


Let  t  3  (x-x') 


•  H  e‘,kVdx  I  •,‘v'  «*.**>*<* 


■  .'-v-'.'*'.'-  v-v v-v *-■.%  v-‘ '-■.yKv v^v- 

vS  v  ‘A*  r  *  s  v-  nk>_-Vr  j* 
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where 


ikhHl1)(kp) 

K(t,x')  =  - — -  ({cosKx  -  cosKx')  +  KxsinKx') 


p  =  [t^  +  h^(cosKx  -  cosKx')^]^  . 


Then, 


1  feikV'  dx'  f  e-,kV  K(,,x')dx 


1  f  e-ik(Van)x'dx-  f  «-,kV  e+fkV'  K(t.x')dx- 


1  f  -1{«-«)Kx'dx,  7  e-1kVK(x,x')dt 


where  k(a  -a„)  *  K(m-n) 
m  n 

and  K(x,x')  is  given  by  Equation  (75).  To  determine  Vm  n,  we  consider  Vm  n  as 

Vm  „  55  Vm  _  +  (V_  n  -  Vm  J  (78) 

171)11  (Df  n  iDjii  m,n 

where  V  is  obtained  using  K(x,  x')  in  the  Equation  (74)  K(x,x')  is  different  fro 
m,n 

K(x,x')  in  that  p  is  set  equal  to  x.  (i.e)  the  cosine  terms  in  p  are  neglected. 

This  manipulation  will  simplify  the  final  computation  of  V  . 

m,n 


0-1  L-i‘Cm-n)Kx'dx-  [  e-1kVK(,,x-)d, 

m,n  a 


.  1  fe-,aKx'x-  [  e-1bTK(x.x')dT 


1  f  ibx  . 

I  j  e  d, 

-CO  0 


f  e“iaKx'K(x,x')dx' 


[• 


K(t,x')  = 


[c(x'+  t)  -  c(x')  -  TC'(x')] 


(80) 


where  a  a  (m-n)  and  b  *  -ka 


m 


Integrating  with  respect  to  x 

ik  H1 


A  kH,(1)(k|T|) 


„  1  ,•  K  -•*  ni  T  ; 

Vn  4|elbTdTl  T  Hrr—  [^(X"T)  - 


C(X-)  -  TC'(x-)]  e"iaKx"  dx- 


(81) 


Let 


1 

A 


c(x')e-iaKx'dx- 


(82) 


Then  from  the  periodicity  of  the  terms  in  the  integrand,  we  have 


Ik  ?  -ibr  r_  /_iaKx 


m 


*n  "  T  J  e 


FT 


Cca(eiaNT-l-1aKT)]  dx 


(83) 


This  integral  is  in  the  form  of  Weber-Schasheitlin  type  (Watson  1944)  and  it 
gives 


Vn  *  V  [(-tr^T7?  * / 1(2  -  <»  +  aK>2  - /  k2-b2J 


(84) 


Substituting  for  a  and  b,  we  have 

(n-m)Ka 


V  *  ic  [■ 
m,n  sm-n  L 


m 


'm 


-  kYm  *  kYn^ 


(85) 


In  the  case  of  a  pure  sinusoid,  of  the  form  hcosKx,  the  fourier  coefficients  are 
h/2  for  a  *  +  1 . 


Vn  ■  1  i  £- 


(n-m)Kc 


m 


'm 


*  k(VYtn>3 


I  m-n I  =  1 


I  m-n I  f  1 


(86) 


6.N.  Watson.  A  Treatise  on  Bessel  Functions,  Cambridge  Univ.  Press 
Cambridge,  England,  1944. 
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The  next  step  is  to  evaluate  (V_  „  -  Vm  )  =  V„ 
r  m,n  m,n'  m,n 


V<im,n  *  1  Je'1l<“"X  dx  I  e+1k“nx'f(x,x')  dx- 


Hi  (ko)  H.(1)(kp) 

f(x,x")  =  ik  p- - r - )  ( (^ (x)  -  s(x"))  -  (x-x'k"(x"))dx- 

2  p 

p  =  [(x-x')2  +  (c(x)  -  c(x'))2]1/2  ( 

P  =  (x-x") 

writing  the  integral  over  the  infinite  limits  as  a  sum 


n  *  j  je”lkamx  dx  l  e1kanx  f(x,x"  +  qA)  dx' 


qa-ooQ 


e‘ikomx  dx 


A 

.  00 

J  e1 kanX"  (  l  f(x,x"  +  qA)  dx") 


A  A 

1  f  -ika_x-imkx  .  f 


A  e  0 


■  eika0x'  +  iknx'(  j  f(XiX,  +  qn) 

« 


=  1  f  e"imkx  dx  |  e1knx'f(x,x')dx' 


where 


f(x,x")  =  e'ikcto(x'x^  £f(x,x"  +  qA).  (90) 

q 

We  can  easily  evaluate  f(x,x")  and  determine  its  fourier  coefficients  corres¬ 
ponding  to  m  and  n.  We  will  show  that,  as  the  point  x-*x"  and  q  =  0,  f(x,x") 
is  bounded.  To  show  this,  we  take  a  sinusoidal  surface  represented  by  c(x)  = 
hcosKx 

H1(1)(kp)  H,(1)(kp) 

f(x,x")  =  (- - - - ; - )  [(h  cosKx  -  h  cosKx")  +  h(x-x")KsinKx"] 

p 

writing  x  =  x"  +  t  as  before 


wwmp 


-  Cm  ** ^  a  •"*  -"o  ’*■  ^ *  ■  .  'C-  ■**,  ■*,  ’ 


w; 
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iHC1>(kp)  ht(1)0<p) 

f(x,x')  ■  (-- - - - ; —  [h  cosK  (x'  +  t)  -  h  cosKx'  +  hKtsinKx']  (92) 

p  p 


as  x-*»x" ,  t-*0 


h[cos(for  +  Kt)  -  cosKx'  +  KtsInKx'] 

*  h[cosKx'*cosKt  -  sinKx'sinKt  -  cosKx'  +  KtsinKx'] 


For  small  values  of  t,  sinKt  »  Kt  and  cosKt  *  (1  -  Therefore,  the  right- 

2  2 

hand  side  of  Equation  93  now  becomes  [-cosKx'(hK' t  /2)].  By  a  similar  argument. 


we  write  for  small  values  of  t 


t(l  +  h2K2sin2Kx')1/2 


Also  for  small  argument 


Hi(1><*>  -  -  3$ 


Substituting  all  these  in  the  expression  for  f(x,x')  we  obtain 


h3K*  cosKx*  sin2  Kx' 


f(x.x-)  .  5ln  £ 

as  x-*or  iwk(l  +  h^K^sIn^Kx') 
q->o 


f(x,x')  can  therefore  be  evaluated  without  difficulty  and  a  2D  FFT  is  performed 


to  obtain  the  elements  of  the  matrix  T 


From  Equation  71 ,  we  have 


[I  -  7  -  7d]$  -  zi 


The  expression  for  Vm  n  has  a  singularity  at  ym  ■  0  and  therefore  we  multiply 
each  equation  by  the  corresponding  y 


[y  -  y7  -  y7d]ij)  ■  2‘Ytp 


where  y  is  a  diagonal  matrix  containing  ym  along  the  diagonal. 


The  method  of  Hoi  ford  can  also  be  used  to  calculate  the  scattering  strength 


due  to  a  beam  of  finite  bandwidth  by  the  following  scheme. 


5S: 


The  incident  plane  wave  is  split  into  its  fourier  components  given  by 


B<“o)  *  3T  {  P1nc‘x-k>e'1k“ond« 

•  00 

The  resultant  scattered  field  is  computed  from 

oo 

'Em  *  \  B(“o)  'W*-z*o>  dao 

•OO 

where  pscat(x,2,ao^  1S  the  scatter  field  due  to  unit  amplitude  plane  wave 
incident  upon  the  surface  with  a  direction  cosine  aQ. 


(99) 


(100) 


III.  Results 

A  possible  model  for  study  of  the  urban  areas  with  a  SAR  would  be  to  model 
the  surface  as  a  rectangular  periodic  surface  with  a  period  of  30m  and  a  height 
of  about  15m.  For  an  incident  wave  length  of  above  .235m,  the  number  of  scattering 
angles  will  be  of  the  order  of  300  and  a  corresponding  300  x  300  matrix  for 
obtaining  the  scattering  intensities  at  these  angles.  Because  of  such  a  large 
matrix  inversion,  a  smaller  period  for  the  surface  was  picked.  It  is  felt  that 
this  scaled  down  model  will  also  exhibit  the  principal  characteristics  of  the 
original  model  we  are  interested  in.  The  computation  involves  estimating  the 
scatter  intensities  at  discrete  angles  given  by  the  grating  equation  using  the 
exact  method  due  to  Holford.  This  was  done  at  different  incident  angles  ranging 
from  5°  to  25°,  with  the  corresponding  back  scatter  directions  given  by  -5°  to  -25° 
respectively.  Higher  angles  were  not  considered  as  bulk  of  the  energy  was 
scattered  in  the  forward  direction  at  these  angles.  Also  the  look  angles  for  SAR 
is  in  the  region  20°  +  3°. 

The  issues  in  the  numerical  implementation  of  the  scheme  are 
(a)  The  infinite  sum  in  Equation  (2.89)  can  be  approximated  by  a  finite 
number  of  terms  and  a  length  of  10A  was  sufficient  to  obtain  con- 
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vergent  results 

(b)  Nyquist  criteria  was  used  for  sampling  while  performing  the  2-D  FFT. 

(c)  A  rectangular  window  for  the  FFT  did  not  perform  well,  while  Hanning 
window  significantly  improved  the  results. 

(d)  Conservation  of  energy  was  used  as  a  test  to  check  the  validity  of  the 
results  and  except  for  one  case,  the  energy  conservation  was  not 
violated.  In  that  case  (5°  angle),  the  reflection  coefficient  was 
slightly  greater  than  one  at  scattering  angles  close  to  grazing 

and  less  than  one  at  other  angles. 

Figures  1  to  4,  the  reflected  energy  is  given  as  a  function  of  the  incident 
angle  (e)  for  various  scattered  angles.  The  back  scattered  angle  is  -e.  As  the 
incident  angle  increases  the  back  scattered  intensity  decreases  significantly. 
Figure  (5-7)  shows  the  back  scattered  energy  as  a  function  of  the  incident  angle 
for  Kh  *  5,10,15  and  A/X  *  10  and  Figure  8  for  Kh  «  TO  and  A/x  »  15.  The  back 
scattered  energy  decreases  rapidly  to  insignificant  levels  beyond  25°.  However 
it  seems  to  be  having  a  minimum  between  0  and  25°.  This  may  be  a  possible 
explanation  for  the  observed  variation  in  the  intensity  of  the  two  SAR  images  in 
the  LA  area. 

For  the  same  Kh  *  10,  as  A/x  increases,  (Figures  6  and  8)  the  reflected 
energy  decreases  significantly.  This  indicates  that  the  fundamental  component 
of  the  spatial  structure  contributes  a  significant  amount  of  energy. 
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IV.  Conclusions  and  Recommendations 

The  numerical  results  indicate  that  there  is  a  definite  variation  in  the 
back  scattered  intensity  as  a  function  of  the  aspect  angle.  However  in  order 
for  us  to  do  a  more  realistic  model,  we  would  need  the  facilities  of  a  large 
computer  like  the  CRAY,  due  to  the  computer  intensive  nature  of  the  calculation. 
In  any  future  work,  a  request  for  such  facilities  has  to  be  made,  and  only  then 
can  one  compare  the  actual  observations  from  the  SeaSat  images  with  the  model 
calculations. 

Modelling  of  the  back  scatter  as  being  from  a  series  of  regularly  spaced 
flat  rectangular  radiators  parallel  to  the  earth  surface  was  not  pursued  as  it 
was  felt  that  significant  errors  would  result  from  neglecting  the  height  of  the 
buildings.  To  include  this  effect,  one  would  have  to  use  some  modified  aspects 
of  ray  tracing  and  GTD  corrections  to  rays.  This  represents  another  alternate 
approach  that  needs  to  be  investigated  and  numerically  implemented. 

In  order  for  us  to  get  a  more  quantitative  feel  for  the  back  scattered 
radiation,  it  would  be  useful  to  perform  experiments  on  scaled  down  models  using 
laser  beams  and  optical  frequencies.  Measurements  of  optical  backscatter  can  be 
made  by  modelling  the  city  as  a  grating  or  a  crossed  grating  and  simulating  the 
satellite  radar  with  coherent  optical  radiation.  The  ratio  of  the  optical 
wavelength  to  the  length  of  the  gratings  periodicity  will  be  equal  to  the  ratio 
of  the  radar's  wavelength  and  the  dimensions  of  a  city  block. 
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